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Thermal emission in the prompt phase of gamma-ray bursts 



F. Ryde 

Stockholm Observatory, AlbaNova, SE-106 91 Stockholm, Sweden 

I discuss the interpretation of the prompt phase in gamma-ray bursts as being dominated by quasi- 
thermal emission, rather than by synchrotron emission. Such an interpretation gives a more natural 
explanation of (i) the observed variety of spectral shape and spectral evolution, (ii) the observed 
narrowness of the distribution of peak energies, as well as (iii) the observed correlations between 
peak energy and luminosity. However, the physical setting that could produce such a scenario is not 
yet clear. 
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I. INTRODUCTION 

The prompt emission of gamma-ray bursts (GRBs), 
radiating mainly as gamma-rays and X-rays, has de- 
fied any simple explanation, despite the presence of 
a rich observational material and great theoretical ef- 
forts. This is in contrast to the afterglow emission, 
in many cases detected all the way from X-rays to 
radio wavelengths, which is successfully described by 
synchrotron emission from a forward shock moving 
at great speed into the surrounding medium. Syn- 
chrotron emission is also a natural candidate for the 
prompt emission since it arises naturally in an ultra- 
relativistic outflow in which the kinetic energy is dissi- 
pated through, for instance shocks or magnetic recon- 
nections, and shared between the magnetic field and 
particles. However, there are several observational 
facts that contradict such a simple picture, most im- 
portantly the existence of spectra which are much too 
hard, see for instance the spectra from GRB930214 
in Figure 1 and GRB960530 in Figure 2. In many 
cases such spectra are fitted well by a thermal emission 
function 0,1113. Furthermore, Ryde 0] showed that 
spectra from more typical bursts, that is, bursts hav- 
ing spectra which are consistent with the synchrotron 
model, can indeed be fitted with a hybrid model which 
is dominated by a thermal component, but that is 
overlayed with a non-thermal emission component as 
well. In many cases such a model gives a statisti- 
cally better fit. Such an example is given in Figure 3. 
Even though bursts appear to have a variety, some- 
times complex, spectral evolutions, the behavior of 
the two separate components is remarkably similar for 
all bursts, with the temperature describing a broken 
power-law in time. The non-thermal component is, 
in most cases, consistent with emission from a popu- 
lation of fast cooling electrons emitting optically-thin 
synchrotron emission or non-thermal Compton radia- 
tion, giving a power-law slope of the photon spectrum 
of s — —1.5. However, in the case of GRB960530, 
shown in Figure 2, s is closer to —2/3, which is ex- 
pected for slow cooling 

It is very important to note that it is the instan- 
taneous spectra of GRBs that most closely should re- 



veal the radiation mechanism. This is because of the 
strong spectral evolution that normally occurs dur- 
ing a burst and that will make the spectral shape of 
the time-integrated spectrum differ from that given 
by the emission process. This is in particular the case 
for complex bursts with several emission peaks. The 
time-integrated spectrum can easily be found from 
the instantaneous spectra and the spectral evolution, 
which was shown analytically by Ryde & Svensson j(| ; 
the spectra always become softer. 



II. SYNCHROTRON MODEL 

It was early recognized that the spectra of gamma- 
ray bursts (GRBs) have a non-thermal character, 
with emission over a broad energy range (e.g. Q)- 
This typically indicates emission from an optically- 
thin source and an initial proposal for GRBs was 
therefore an optically-thin synchrotron model from 
shock-accelerated, relativistic electrons (e.g. H, EJ)- 
The number density of the radiating electrons is as- 
sumed to be typically a power law as a function of the 
electron Lorentz factor "f c above a minimum value, 
7min, with index —p. Such a distribution gives rise 
to a power-law photon spectrum with photon index 
a — —2/3 below a break energy E p oc j^ in and a 
high-energy power-law with index (3 = — (p + l)/2. 
However, as mentioned above, this model has difficul- 
ties in explaining the observed spectra of GRBs, which 
show a great variation in a and j3 (see ^3)- ^ n par- 
ticular, a substantial fraction of them have a > —2/3, 
which is not possible in the model in its simplest form, 
since a = —2/3 is the power-law slope of the fun- 
damental synchrotron function for electrons with an 
isotropic distribution of pitch angles [Tlj . The prob- 
lem becomes even more severe for the case when the 
cooling time of the electrons is shorter than the typi- 
cal dynamic timescale. In the typical setting of GRBs 
having a relativistic outflow with a bulk Lorentz fac- 
tor F ~ 100, the time scales for synchrotron and 
inverse Compton losses are ~ 10 -6 s [13, which 
is much shorter than both the dynamic time scale 
R/2T 2 c ~ 1 s (i?/10 15 cm), and the integration time- 
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FIG. 1: Two time-resolved spectra from GRB930214 (BATSE trigger 2193) from 6 and 40 seconds after the trigger. 
Note that the spectra are fitted well with a Planck function, both in the Rayleigh- Jeans portion of the spectrum with 
a — +1 and in the Wien portion with a fast voidance of flux. The temperature has changed between the measurements. 



scale of the recorded data, typically 64 ms to 1 s. In 
such a case the low- energy pow er-law should be even 
softer, with a = —1.5 |l3l Il4|. now contradicting a 
majority of the observed spectra. Furthermore, the 
observed distribution of a from instantaneous spectra 
is smooth (see 01) an d does not show any indica- 
tion of preferred values, such as —2/3 or —3/2. Other 
variations of the synchrotron or/and inverse Compton 
model have been suggested (see e.g. 0, 0, 0j) to 
account for these hard spectra. 

The peak energy from the above distribution of elec- 
trons is given by E p ^ — 7 m £?j_r. In the external shock 
model 7 m and B±_ are proportional to the bulk Lorentz 
factor, which makes Ep^ oc T 4 , which is a very strong 
dependence, which poses a problem in explaining the 
relative narrowness of the distribution of peak ener- 
gies ^(J , even including the X-ray flashes. For the 
internal shock model the 7 m oc r rc i, which is the rela- 
tive Lorentz factor between the two shells that collide 
and 



E pk oc B±_T. 



(1) 



However, the sharing of the energy between the kinetic 
energy of the electrons and the magnetic fields should 
lead to a larger dispersion. 

A third complication arises in explaining the ob- 
served correlation between peak energy and the lumi- 
nosity which was discussed by Lloyd-Ronning et al. 
[l8| and Amati et al. 0] (see also H3); the peak 
energy is correlated with the isotropically equivalent 
energy given by 



E, 



P k 



Eh, 



0.40±0.05 



where E\ sa is 



1.2 x 10 53 erg 



E 1 = (l- cos9)Ei, 



(2) 



(3) 



where E n is the actual gamma-ray energy emitted and 
9 is the jet opening half-angle of the collimated out- 
flow. 

Equation shows that in the internal shock model 
Epk is proportional to the Lorentz boosted magnetic 
field strength. The total energy density 



U 



(BTf L_ 



(4) 



The typical radius for the internal shocks to occur is 
i? s h ~ ci v r 2 , where t v is the typical variability time 
scale, and thus 



£ pk cx T- 2 L 



1/2.-1 

L 1T 



(5) 



To get a relation similar to that in equation J5J both 
r and t v have to be quite similar for all bursts, which 
is difficult to imagine. Even though there is no di- 
rect way of determining the bulk Lorentz factor, var- 
ious physical models give suggestions on plausible re- 
lations between the luminosity and the bulk Lorentz 
factor. For instance, 23] argued for L k T 2 . Such 
a relation would thus give i? p k oc L^ 1 / 2 , that is, an 
anti-correlation, in contradiction to the observed be- 
havior (see also Additional assumptions are 
needed to explain the positive correlation. Invoking 
Poynting flux and/or pair dominated models the cor- 
relation also becomes positive pRj . 



III. QUASI-THERMAL MODELS 

If the prompt phase is indeed dominated by a ther- 
mal component these three issues become natural con- 
sequences. First, as shown by Ryde (2004, 2005) the 
relative strength and the slope of the non-thermal 
component will determine the value of the low-energy 
power-law index, a, that would be found if the Band 
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FIG. 2: The same spectrum from GRB960530 (#5478; 6 s after the trigger) fitted with (left panel) the Band et al. 
[2^1 model with a = 1.7 ± 1.5 and /3 = —2.4 ± 0.3 and (right panel) the hybrid model 0], with a power-law slope of 
s = -0.62 ± 0.27. Note the obliging of the data points pl l lS ^ . 
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FIG. 3: The same spectrum from GRB911031 (#973; 2.5 s after the trigger) fitted with (left panel) the Band et al. [2eJ 
model with a = —0.86 ± 0.12 and /3 = —2.4 ± 0.3 and (right panel) the hybrid model with a power-law slope of 
s = -1.52 ±0.04. 



et al. [26( function were to be used. If the thermal 
component is strong and/or the non-thermal compo- 
nent is hard, the resulting spectrum will have a hard a 
(see Fig. 2). While if the non-thermal component be- 
comes relatively stronger and/or softer the measured 
a- value would be softer (see Fig. 3). The observed dis- 
tribution of a-values is therefore consistent with this 
picture and in particular the spectra beyond the " line- 
of-death" |23| are not conspicuous. A strong spectral 
evolution, for instance, in a large change in the mea- 
sured value of a, is also easily explained. 

Second, the peak of the spectrum is now determined 
of kT and is less sensitive to the bulk Lorentz fac- 
tor. In fact, if the photosphere occurs during the 
acceleration phase it is practically independent of V. 



Rees & Meszaros 28] suggest a model where the pho- 
tospheric emission can become enhanced by dissipa- 
tive effects below the photosphere (magnetic recon- 
nections, shocks) and subsequent Comptonization, see 
also j2^| . Typical values for the peak energy would be 
hundreds of keV. 

Third, the correlation in equation Q has a nat- 
ural explanation since for a thermal emitter the lu- 
minosity and the temperature are correlated. For in- 
stance, equating _E p k with the energy density, one gets 
E pk cx TfeT oc TU 1 / 4 oc TiL/TR 2 ) 1 / 4 . Using that the 
pair photosphere occurs at i? p h oc LT~ 3 (see Rees & 
Meszaros |2S|) then E pk oc L~°' 8 . In the last step, 
L oc r 2 was again assumed. Similarly if the pho- 
tosphere is emitted during the acceleration phase its 
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temperature will be constant in the observer frame 
since the comoving cooling by adiabatic expansion is 
compensated for by the increase in T. 

/ r \i/4 

TkT oc i-^j oc L 1 /4 jR -i/2_ (6) 

Here, Rq is the radius at which the linear acceleration 
starts. Assuming, for instance, that i?o oc L^ 1 then 
again oc L 75 (see further Rees & Meszaros 28]). 
Finally for extremely photon starved plasmas E p ^ oc 
L/Ny oc L. 

IV. CONCLUSIONS 

The radiative efficiency of the thermal emission can 
in plausible scenarios be radically increased by, for in- 
stance, dissipation processes below the photosphere 



[1] Ghirlanda, G., Celotti, A., & Ghisellini, G. 2003, 

A&A, 406, 879 
[2] Kancko, Y., Preece, R. D., & Briggs, M. S. 2003, AAS, 

203, 8004 
[3] Ryde, F. 2004, ApJ, 614, 827 

[4] Ryde, F. 2005, ApJ Letters, in press 

( astro-ph/0504450 ) 
[5] Ryde, F., & Battelino, M. 2005, submitted 
[6] Ryde, F. & Svensson, R. 1999, ApJ, 512, 693 
[7] Fishman, G.J., & Meegan, C.A. 1995, ARAA, 33, 415 
[8] Katz, J. I. 1994, ApJ, 432, L107 
[9] Tavani, M. 1996, ApJ, 466, 768 
[10] Preece, R. D., Briggs, M. S., Mallozzi, R. S., Pendle- 
ton, G. N., Paciesas, W. S., & Band, D. L. 2000, 
ApJSS, 126, 19 
[11] Pacholczyk, A. G. 1970, Radio Astrophysics (San 

Francisco: W. H. Freeman and Co.) 
[12] Ghisellini, G., Lazzati, D., Celotti, A., & Rees, M. J. 

2000, MNRAS, 316, L45 
[13] Bussard, R.W., ApJ, 284, 357 

[14] Ghisellini, G., Celotti, A., Lazzati, D. 2000, MNRAS 
313, LI 

[15] Baring, M.G., & Braby, ML. 2004, ApJ, 613, 460 
[16] Lloyd-Ronning, N. & Petrosian, V. 2000, ApJ, 543, 



psl 13 ■ These processes would naturally produce 
large amount of electron-positron pairs with modest 
Lorentz factors, which would Compton up-scatter the 
thermal radiation. The observed peak would then be 
this Comptonized peak. The non-thermal emission 
seen in the spectra, could be due to synchrotron emis- 
sion or inverse Compton emission from dissipation re- 
gions outside the photosphere. Ryde 3 showed that 
the energy flux in the thermal and the non-thermal 
components are correlated which might indicate the 
latter. 

In summary, thermal emission could indeed domi- 
nate over non-thermal emission in standard settings of 
a GRB jet. In such scenarios a correlation between the 
peak energy and the luminosity naturally arises, the 
details somewhat depending on the dissipation pro- 
cesses. In addition, the dispersion in Ep^ would be 
smaller and the observed spectral shapes and spectral 
evolution get natural explanations. 



722 

[17] Stern, B. & Poutanen, J. 2004, MNRAS, in press 
[18] Lloyd-Ronning, N. M., Petrosian, V., & Mallozzi, R. 

S. 2000, ApJ, 534, 227 
[19] Amati et al. 2002, A&A, 390, 81 

[20] Ghirlanda, G., Ghisellini, G., & Lazzati, D. 2004, 
ApJ, 616, 331 

[21] Fenimore, E. E., Klebesadel, R. W., & Laros, J. G. 
1983, in Gamma-Ray Astronomy in Perspective of Fu- 
ture Space Experiments (New York: Pergamon), 201 
[22] Bromm, V. & Schaefer, B. 1999, ApJ, 520, 661 
[23] Kobayashi, S., Ryde, F., & MacFadyen, A. 2002, ApJ, 
577, 302 

[24] Ramirez-Ruiz, E., & Lloyd-Ronning, N.M. 2002, 

New A, 7, 197 
[25] Zhang, B. & Meszaros, P. 2002, ApJ, 581, 1236 
[26] Band, D., et al. 1993, ApJ, 413, 281 
[27] Preece, R.D., Briggs, M.S., Mallozzi, R.S., Pendleton, 

G.N., Paciesas, W.S., & Band, D.L. 1998, ApJ, 506, 

L23 

[28] Rees, M. J., & Meszaros, P. 2005, submitted 
[29] Pe'er, A., & Waxman, E. 2004, ApJ, 613, 448 



0107 



